Strongly reduced gap in the zigzag spin chain with a ferromagnetic interchain coupling 
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We study a spin 1/2 Heisenberg zigzag spin chain model 
near decoupled two chains. Taking into account a symmetry 
breaking perturbation, we discuss the existence of an energy 
gap in the ferromagnetic interchain coupling as well as the an- 
tifcrromagnetic one. In the ferromagnetic model, a marginally 
relevant fixed line reduces the gap strongly, so that the cor- 
relation length becomes an astronomical length scale even in 
order 1 coupling. This result agrees with density matrix renor- 
malization group results. 



I. INTRODUCTION 

It is interesting to study the effects of a frustration in 
a quantum system. The zigzag spin chain is one of the 
simplest quantum spin models with a frustration. An ex- 
periment on the compound SrCuC>2 reports a remarkable 
suppression of the three dimensional ordering tempera- 
ture due to a strong quantum fluctuation enhanced by 
the frustration!!]. Here we study the effect of frustration 
in the zigzag chain with the following Hamiltonian 

H = [MSi ■ Si+i + Ti • T,+i) + JiS, • (T< + T i+1 )} . 

i 

(i) 

The operators Si and Ti are s — 1/2 spins. Here we 
treat an interchain coupling J\ as a perturbation to the 
two decoupled antiferromagnetic (Jg > 0) spin chains. 

It is well-known that the zigzag chain model with an- 
tiferromagnetic region J\ > has a dimer phase and 
gapless phased The region of the dimer phased (0 < 
J\ < J c ~ 4.15 Ja) contains a Majumdar Ghosh point 
J\ = 2J 2 where the dimerized ground state is obtained 
exactlyu. The point J\ = is a critical point of the 
two decoupled antiferromagnetic chains. There is an- 
other critical point J\ = — AJ2, where the level crossing 
between singlet and fully polarized ferromagnetic ground 
states occurs. There is an exact solutioro of a fully po- 
larized ferromagnetic ground state in J\ < — 4J 2 , and for 
a doubly degenerated ground state at the critical point 
J\ = — 4J2. A numerical analysis in a region J\ > — AJ2 
indicates a complicated size dependence of the ground 
state energyBQ. nn 

It has been long believed that this model is gaplesafm 
for a small ferromagnetic region J\ < 0. In several papers 



however, one loop renormalization group (RG) shows an 
instability of the critical point J\ = both in ferromag- 
netic (Ji < 0) and antiferromagnetic (Ji > 0) region 
due to a Lorentz symmetry breaking perturbationOE-j. 
This fact is a puzzling because the expected energy gap 
produced by the unstable flow has never been observed 
in the ferromagnetic region either numerically or experi- 
mentally. If it were gapless, there should be a new stable 
critical point missed in the one loop approximation. In 
this paper, we clarify a natural mechanism to solve this 
problem. We conclude that the actual energy gap is fi- 
nite but very tiny in an extended region of the coupling 
constant space. We can understand unstable the RG flow 
and numerical analysis consistently. The energy gap is 
too small to observe in any numerical method. 

A stable large scale reduction is an important and diffi- 
cult fundamental problem in theoretical physics. In par- 
ticle physics, any unified theory has this kind of hierarchy 
problem. Generally speaking, a unified theory has only 
one large energy scale, and also should explain the gen- 
eration of a small energy scale to describe the low energy 
physics. These two requirements makes unified theories 
quite unnatural, since we need a fine tuning of the cou- 
pling constants for large scale reduction. We believe that 
nature do not chose special value of the coupling con- 
stants. From the view point of hierarchy problem in the 
theoretical physics, this zigzag chain model can be an 
interesting example in statistical physics. We conclude 
that the very tiny gap is always obtained in an extended 
region of the ferromagnetic interchain coupling unlike the 
antiferromagnetic coupling. 

This paper is organized as follows. In section II we 
construct an effective field theory for the zigzag chain 
and calculate the beta functions. In section III, we study 
the RG flow and discuss an energy gap reduced strongly 
in the ferromagnetic interchain coupling region. In sec- 
tion IV, we calculate the energy spectrum to describe the 
results of numerical calculation. In section V, we show 
our numerical analysis. 



II. EFFECTIVE FIELD THEORY 

The unperturbed theory is two decoupled Heisenberg 
antiferromagnetic chains whose-effective theory is two de- 
coupled SU{2)\ WZW modelsEl. We introduce two free 
bosons v? a (z, z) = ifa(z) + (fia(z) (a = 1, 2) with 



So 



2tt 



d 2 z [dipidipi + CV2CV2) 



(2) 



1 



with two point functions (ip a (z)(pb(0)) = — <5 a hlogz, 
((p a (z)(pb(0)) = -Sab log 2, for two WZW models. The 
spin operator is written in terms of two bosons 



2tt S, = Ji + Ji + (-lYMtr[{ gi + g\)—l 



(3) 



2vr Tj = J 2 + J 2 + (-iyMtr[(g 2 +gl)-^}, (4) 

where M is a nonuniversal real constant. Here, the SU(2) 
currents and primaries are written in two bosons and 
Klein factors 

J+(z)=r,air l aie i ^ z \ J z a (z) = ^d<p a (z), (5) 

.9a tt = VaTVa^+^V^, g aU = Varfkie^'-^y^, 
a t = n .p-^a-VaJ/Va „ _ „ ~ -i(<p a +v>a)/V2 

i/aj.y — '/oJ,'/aT e ' ya-ll — '/o|'/ai e ! 

with a = 1,2. The Klein factors obey anticommutation 
relation {rj aa , r)b@} = 25 a b&af), to satisfy the following cor- 
rect operator product expansion for the SU(2) symmetry 



Ja( z )9baf}{w,w) ~ Sab ^ 2 ( t a fc g a ) a/? (w, w), 

z — w 

S /2 

The interaction operators 



d 2 z 

Z7T » — ' 

1=1 



(7) 



can be represented in terms of two SU(2) currents and 
primary fields of the WZW model 

4>i(z, 2) = Ji(z) ■ Ji(z) + 3 2 {z) • J 2 (z), 
4> 2 (z,2) = 3\(z) ■ 3 2 (z) + 3 2 (z) • 3 i(z), 
(fo(z, 2) = tr[gii(d - B)g 2 ], 
<t>i(z, 2) = trgii(d - d)trg 2 , 

fa{z,z) = 3t(z) ■ 3 2 {z) + 3t(2) ■ 3 2 {2), (8) 

where A 2 = Ji/ir, A 3 = -JiM 2 /tt, A 4 = JiM 2 /(2ir), 
A5 = Ji/(27r) and Ai is a certain negative constant. 
The initial coupling constant Ai is roughly estimated 
as the-prder 1/10 in numerical analysis for single linear 
chainl!3. Note that the the perturbations ^3(2,2) and 
4>a(z,z) include operators with the conformal dimension 
(§,5) ancLLi^g), which break the Lorentz and the parity 
symmetryliJ'EJ. The original spin model has SU(2) sym- 
metry, translational symmetry and a permutation with a 
translation of the one chain 



-i+l> 



(9) 



which forbids interactions with dimension x < 2. If this 
symmetry is broken, the symmetry breaking can be mea- 
sured by the following order parameter 



Si • (T< - T i+1 ) ~ tr[( 9l + <?i)y] -tr[(g 2 

(10) 

This order parameter also changes sign under a trans- 
lation of one spin Tj — > Tj+i. The operator product 
expansion 

&(2,2)&(Q > 0)~5^^&(0,0), 
k |z| 

gives one loop renormalization group beta function 
.dX k 



-l- 



dl 



9 CijkW 



which has 



he following practical form 
Ax — A3A4 — A4, 



~dT 

dA, 

dl 
dh 

dl 
e?A 4 

~dT 

dX 5 

dl 



— A 2 . + A3A4 + A3, 
1 3 

3 1 
= A1A3 + 2^!^ 4 — 2^ 2/ ^ 4 ' 

= 7^X3X4- 



(11) 



Note that the beta functions have a symmetry A3 — ► —A3, 
A4 — * — A4, which corresponds to a translation on one 
chain T, — > T i+1 in the original spin model. 



III. RENORMALIZATION GROUP FLOW 

The critical point J\ = divides ferromagnetic and an- 
tiferromagnetic dimer phases. Both regions have energy 
gap. The renormalization flow diverges eventually along 
a stable direction (Ai, X 2 , A3, A 4 ) = (0,1,-^,0) for the 
ferromagnetic coupling J\ < 0, and along another stable 
direction (Ai, X 2 , A3, A4) = (0,1, ^=,0) for the antiferro- 
magnetic coupling J\ > 0. The current-current coupling 
A 2 grows with positive value both in ferromagnetic and 
antiferromagnetic region. In the effective field theory, 
the ferromagnetic model does not differ from the anti- 
ferromagnetic model transformed by Tj — > Ti+i, which 
changes the sign of A3 and A4. This fact suggests that 
both ferromagnetic and antiferromagnetic models have 
the same dimerization pattern. The correlation length £ 
behaves as 



£ ~ aexpcj J\\ 



(12) 



with v = 2/3 for J\ ~ and the lattice spacing a. This 
result iS|-obtaincd in another RG method for the RG 
equatioril3. In the antiferromagnetic region, this scal- 
ing formula of the correlation length holds for relatively 



2 



large value in < J\ < J c ~ 4.15J2, even though the 
formula is obtained merely by one loop approximation. 
We demonstrate it by numerical results from density ma- 
trix renormalization group calculation(DMRG)tj in Fig. 
[l]. The gap scaling formula can fit the numerical result 
quite well. This agreement even in relatively strong cou- 
pling region is not so surprising. There are many systems 
with Kosterlitz Thouless type phase transition, where the 
gap scaling formula calculated in the weak coupling re- 
gion holds even in the strong coupling region. In the 
ferromagnetic region, however, this scaling formula does 
not hold for | J\ | ~ J 2 ■ The energy gap in this region 
cannot be found in any numerical analysis. There should 
be some special reason. 
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FIG. 1. The energy gap A vs. the antiferromagnetic inter- 
chain coupling Ji > is depicted. The solid line is drawn by 
the scaling formula A ~ exp(— c| Ji \~ 2 ^ 3 ) for small Ji in the 
infinite order phase transition, using Eq.(p|) with £ ~ 1/A. 
The gap A calculated by the DMRG is obtained from Ref. [ 

§■ 

We clarify in the following why the gap scaling formula 
in ferromagnetic region differs from the antiferromagnetic 
region. The beta functions for Ai, A2, A3, A4 have simul- 
taneous zero on a line 

Ai = A 2 = 0, A 3 + A 4 = 0. (13) 

This line behaves like a fixed line for a four dimensional 
coupling constant space (Ai,---,A4), since A5 does not 
enter their beta functions. Near this line, we define the 
deviation of the running coupling constants from the line 
by A3 = A + SX3, A4 = —A + 6X4. The linearized beta 
functions are 

4(<*Ai-SA 2 )~0, 
al 

Ij^i + 6X 2 ) ~ 2X(SX 3 + 6X4), 
l^(5X 3 + 6X 4 )~-X(6X 1 -dX 2 ), 



1^(6X3-6X4)^0, (14) 

The eigenvalues of the scaling matrix all vanish on that 
line. This marginal property of the fixed line yields a 
remarkable phenomenon. The linearized flow near the 
fixed line can be integrated as follows 

6Xi(l) - 5X 2 (l) ~ A, 

6X x (l) + 5X 2 (l) ~C + 2XB\nl 

6X 3 (l) + 5X 4 (l) ~ B- XAlnl 

6X3(1) - 6X4(1) ~ D, (15) 

where A, B, C and D are integration constants deter- 
mined by initial coupling constants which are nonuni- 
versal. If the initial coupling constants lie near this fixed 
line, A + C is negative, A - C ~ -2Jx/-k, B + D = 
- JiM 2 /tt - A and B - D = J 1 M 2 /(2tt) - A. We can 
study the nature of the RG flow numerically together 
with an analytic argument based on the behavior of the 
fixed line. Although the line is unstable except the case 
A = B = 0, there is an extended region where the run- 
ning couplings 6Xi(l) + 6X 2 (l) and 6X3(^+6X4(1) flow to- 
ward 0. In this case other couplings are renormalized log- 
arithmically, and the running coupling constants spend 
a long time (long length scale) near the fixed line (]13]). 

In the ferromagnetic region with small we have 
A < 0, B > 0, C < 0, A > 0. In this case, 6X3(1) + 6X A (l) 
does not flow toward 0, and then the flow is free from 
the fixed line as well as in the antiferromagnetic region. 
In this region, we consider the gap scaling formula ( |l2] ) 
holds and the gap is still small enough. At J\ ~ — O.IJ2 
however, the constant A changes sign and the running 
coupling constants start to flow toward 0. The corre- 
lation length grows again by the effect of the fixed line 
(|l3|). The correlation length can be estimated approx- 
imately from the length scale / where the running cou- 
pling constant Xi(l) diverges. The numerical solution of 
the RG equation shows a dependence of the correlation 
length on the non-universal constant M. Typically, the 
minimal correlation length in the ferromagnetic region 
becomes an astronomical length scale more than 10 36 a 
for J \ ~ — O.2J2 and M = 1. The interaction between 
two chains may change the non-universal quantities Ai 
and M 2 . However, the qualitative nature of the flow is 
unchanged for M 2 < 2 and J x < 0. If M 2 > 2.5, the gap 
may be observed by numerical analysis. We estimate the 
gap of the zigzag chain by numerically evaluating the 
renormalization group equations. The data In A vs. J\ 
is plotted in Fig. ^| setting the nonuniversal quantities 
AI = 1 and Ai = — 0.24/-7T as typical values. Also in this 
case we could observe asymmetric property between the 
ferromagnetic and antiferromagnetic coupling J\. The 
theoretical fitting function with v = 2/3 can fit the data 
in the antiferromagnetic side J\ > better than those in 
the ferromagnetic side J\ < 0. 
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Fig. 2, Itoi et al, PRB 
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FIG. 2. The logarithm of gap In A vs. small Ji estimated 
by the scale where the solution of the RG equation has the 
singularity. The fitting line is obtained by fitting formula 
A = a exp(— c/| Ji | 2//3 ) with data close to zero. Gaps obtained 
from RG equations in between 0.1 > Ji > —0.1 is used in 
fitting. 

This fact implies a quite unusual phenomenon. The 
infinite system differs essentially from a macroscopic sys- 
tem with the finite size. A system with a macroscopic 
finite size is described in a massless theory on the fixed 
line, while an infinite system is described in a massive 
theory. We can employ the massless field theory for a 
macroscopic system available in an ordinary condensed 
matter experiment, since the correlation length becomes 
an astronomical length scale. We can understand the 
quite slow convergence in numerical methods in this re- 
gion. Although th e bet a functions depend on the renor- 
malization schemet2lliil, the existence of the marginally 
relevant fixed line is scheme independent and the strong 
gap reduction in an extended region occurs universally. 
On the other hand, in the antiferromagnetic interchain 
coupling, the flow is free from the fixed line and the cor- 
relation length becomes £ ~ 7a for J\ ~ J-2- This result 
is consistent with the numerical analysis in Ref.[ ^| as 
depicted in Fig. |l|. 



IV. LOW ENERGY LEVELS 

Next we consider the effect of the chiral operators 
4>5(z) = Ji(-z) • J2(z) and <fe(.z) = Ji(z) • Ja(^) with con- 
formal dimension (2, 0) and (0, 2) respectively. In an ordi- 
nary macroscopic system size, we can employ the mass- 
less theory as a good effective theory. A real fermion 
representation is useful to see the effect of the chiral 
operatoioliSrt 



€i(z) + i&(z) = W-77+exp 



ipi(z) + <p 2 (z) 



£3(2) + i&(z) = \ -V- exp 
a 



V2 



v/2 



(16) 



Here r}± are written in a bilinear of rj aa such that the total 
current for the SO(3) generator becomes Ji + J2 = §£x£, 
where £ = (£1, £2, £3). Then the fields £ and £4 corre- 
spond to a triplet and singlet excitation, respectively. In 
this mapping, the chiral perturbation operator has a free 
field representation 

cj) 5 (z) = --dipidip 2 + V11VHV21V21 coB[V2(ipi(z) - <P2(z))] 



(17) 



The chiral perturbation operator ([Tt]) makes the shift 
of the spin wave velocities of triplet and singlet fieldsO 



Vt 



Ji 
4tt' 



3Ji 

47T 



(18) 



These shifts are found in the finite size correction to 
the low energy levels of two triplet boundary operators 
£j£fc, £i£i with j,k — 1,2,3 in open boundary condi- 
tion(OBC) and two triplet operators (Jj^j^k^A, o'il^jl^k^A 
for periodic boundary condition (PBCp. Each gap in the 
leading finite size correction in OBC depends on J\ 



AE. 



AE, 



irv t 7r(u + Ji/Air) 
^ ~T = L ' 

ir(vt + v s ) _ vr(u - Ji/47t) 



(19) 



3 2L L 

In PBC, all belong to the same energy independent of J\ 



AE. 



27r(3ut + v s ) 2nv 
4i ~ ~T~' 



(20) 



There are some logarithmic corrections for a short chain 
mainly through Xi(L) + \2(L) ~ 2/lni much larger than 
the effect of renormalization SX^L) ~6x 10~ 4 lnL. We 
can expect the flat band instability at J\ — —A~kv in 
the fermionization. This instability is a sign of a level 
crossing between the singlet and the ferromagnetic states. 



V. NUMERICAL ANALYSIS 

These arguments are consistent with DMRG analysis 
on the OBC zigzag chain. Here we chose J2 = 1 and 
calculate-, the low energy levels numerically by DMRG 
methocO. The gap is too small and the correlation 
length is too large to check the energy gap directly for 
Ji < in numerical analysis. However, the results ob- 
tained in previous sections can be supported by the low- 
est three energy levels of the system. The ground state 
and first excited state for each chain at J\ = are unique 
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Fig. 4 Itoi et al, PRB 



and of total spin zero and one respectively^ for even 
length chains. We label the ground state for the zigzag 
chain as (0,0) which is the product of the ground states 
of the two chains. The two degenerate first excited states 
of total spin one, (0, 1)' and (0, l) s , are the parity odd 
and even states as the product of the ground state and 
the first excited state from different chains. The parity 
is for the symmetry of S chain and T chain permutation. 

When Jj decreases from zero, the excitation energy 
behaves as gapless for the infinite system in our finite 
size scaling fitting. The energy of (0,1)*, irvt/L de- 
creases, and the one ir(v s + v t )/2L of (0, l) s increases 
for small \Ji\. In Fig. [| we demonstrate the two dif- 
ferent spin velocities obtained by scaling for an arbitrary 
point Ji = —0.5 in the phase with small and negative J%. 
We have kept m = 500 states in DMRG calculation for 
OBC chains, and the biggest truncation error in DMRG 
is 10" 7 . 

Fig. 3 Itoi et al, PRB 
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FIG. 3. The spin velocities for Ji = —0.5. The excitation 
energies calculated by DMRG for the state (0, 1)' and (0, l) 3 
multiplied by L is plotted vs. 1/lnL with L up to 48. Two 
different spin velocity v t and v s are obtained for them respec- 
tively by least square fitting. The system behaves as its gap 
scaling being 1/L. 

The linear increasing of vt and decreasing of v s pre- 
dicted in Eq.(^8|) for J\ ~ can be observed in Fig. [l|b 
the low energy levels for OBC L — 8 chains. In Fig 
DMRG result is exact for such a short length L = 8. We 
have plotted the lowest five energy levels for each S l ° tal 
in the figure vs. the interchain coupling J\. 

In Fig. ||, the low energy levels obtained by exact di- 
agonalization for a PBC L = 6 chain is shown. We can 
check in the figure that the two lowest excitation energies 
are independent of J±, as we expected in Eq. (p?j|) . The 
two lowest excitation energies are of total spin 1. In the 
figure, eight lowest energy levels for each total momen- 
tum and Sl otal of L = 6 (12 sites) PBC chains has been 
plotted vs. interchain coupling J\. 
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FIG. 4. Low excitation energies calculated by DMRG for 
zigzag open chain for length L = 8 with negative Ji. The 
lowest excited state at Ji ~ is (0, 1)' and the second lowest is 
(0, l) s . The states collapse to zero energy when Ji approaches 
—4. A non singlet ground state seems to appear at Ji ~ —2. 
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FIG. 5. Low excitation energies calculated by ex- 
act-diagonalization for zigzag PBC chains with L = 6 (12 
sites). For each negative Ji, eight lowest energy levels for 
each total momentum and sl° tal is plotted. The two lowest 
excited states at Ji ~ are of total spin 1. Their energy 
is degenerate and almost independent of Ji as predicted in 
Eq.©. 

The low energy excitations are described in the crit- 
ical theory near J± — with the shifted spin wave ve- 
locities, therefore we can naturally exclude a different 
critical theory at a new fixed point away from J\ = 
which describes this system. 

In Fig. |] flat band instability is demonstrated. It oc- 
curs at J i = — 4 and the system enters the ferromagnetic 
phase with the maximum spin S tot = L ground i^tete, 
which is consistent with previously obtained resultsafi. A 
non singlet ground state,seems to occur near J\ = —2 as 
pointed by Cabra et alEJ. It is difficult to judge whether 
this partially polarized ground state can survive in the 
infinite length limit. In a macroscopic system size, it 
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seems alive in J\ > —4. 

VI. DISCUSSIONS 

In conclusion, the divergent RG flow by the Lorentz 
symmetry breaking perturbation certainly produces an 
energy gap in the zigzag chain with the ferromagnetic 
interchain coupling as well as in the antiferromagnetic 
model. In an extended region of the ferromagnetic inter- 
chain coupling with order 1 , the correlation length can be 
extremely large due to the effect of the marginally rele- 
vant fixed line. If there is a finite energy gap, it becomes 
too tiny to observe. The quantum fluctuation enhanced 
by the frustration yields the extraordinary reduction of 
the energy gap, which cannot be checked in ordinary 
macroscopic physics. The slow convergence in numerical 
analysis can be understood for this reason. The DMRG 
analysis shows that the massless theory near the two de- 
coupled chains describes well the low energy levels in a 
macroscopic system. 

Here, we comment on an experimental result on the 
zigzag chain compound SrCu02 and the linear chain com- 
pound Sr2CuC>3. Motoyama, Eisaki and Uchida show 
that the temperature dependence of the spin susceptibil- 
ity of the zigzag chaip-,SrCu02 differs from that in the 
linear chain Sr 2 Cu03li3. The drastic decreasing of the 
susceptibility is observed in the linear chain Sr 2 Cu03 in 
low temperature. This phenomenon is not observed in 
the zigzag chain SrCu02, even though the three dimen- 
sional ordering temperature ~ 2 K of the zigzag chain 
SrCuC>2 is less than the temperature ~ 5 K of the linear 
chain Sr2Cu03. We would like to point out the fact that 
the renormalization group flow (pT|) has a certain param- 
eter region of the ferromagnetic interchain coupling in 
which 

x(T ) ~ -1 (1 - Ax(l/T)/4 - A 2 (l/T)/4) (21) 

does not decrease so drastically in the low temperature 
region. It might be possible to understand the difference 
of the low temperature behavior of the susceptibility be- 
tween the linear chain and the zigzag chain in terms of 
RG flow. 

The universal scale reduction in the zigzag chain is in- 
teresting as a rare example from the view point of the 
hierarchy problem in the theoretical physics. A string 
theory as an ultimate unified theory of every thing in 
particle physics needs a natural explanation of the strong 
suppression of the energy scale. Such a unified theory 
should explain the hierarchy of masses of the elementary 
excitations as well as the hierarchy of all interactions. We 
should obtain the mass of the light particle as a reduction 
from the Planck mass ~ 10 28 eV, for example, the neu- 
trino mass ~ 1 eV, the electron mass ~ 10 6 eV, the muon 
mass ~ 10 8 eV, the Z boson mass ~ 10 11 eV. To under- 
stand this mass hierarchy problem in a unified theory, 



we have to answer the following two questions. First, 
why do several energy scales appear depending on the 
species of the fundamental excitations ? Second, how is 
the energy scale reduced with very large rate? Normally, 
the small mass of fcrmionic excitations is considered as 
the result of a small chiral symmetry breaking. For the. 
small mass of bosonic excitations, the super symmetryEJj 
or composite approaches may work. For the too large 
reduction rate, however, we have to tune the coupling 
constants to the corresponding small region. This is pos- 
sible, but unnatural. These two questions seem too dif- 
ficult to resolve. In condensed matter physics, however, 
we have interesting examples to consider these problems. 
For the first question, we can refer the charge-spin separa- 
tion in the Tomonaga-Luttinger liquid. Even though this 
model has the explicit symmetry breaking, the two fun- 
damental excitations have different energy scales in the 
long distance physics. For the second question, a stable 
large scale reduction occurs in the ferromagnetic zigzag 
chain model without any fine-tuning of the coupling con- 
stant J\. In the almost models, as in the antiferromag- 
netic zigzag chain model, however, the fine tuning of the 
coupling constant to the small region realizes the large 
scale reduction, even if the spontaneous chiral symmetry 
breaking works there. The experience on these examples 
may be useful to find a solution of the mass hierarchy 
problem. 

This model should be studied still in more accurate 
approaches than the one loop renormalization group. 
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